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The Deep Space Network (DSN) routinely measures Earth orientation in sup- 
port of spacecraft tracking and navigation using very long- baseline interferometry 
(VLBI) with the deep-space tracking antennas. The variability of the most unpre- 
dictable Earth-orientation component , Universal Time 1 (UT1 ), is a major factor in 
determining the frequency with which the DSN measurements must be made. The 
installation of advanced Global Positioning System (GPS) receivers at the DSN 
sites and elsewhere may soon permit routine measurements of UT1 variation with 
significantly less dependence on the deep-space tracking antennas than is currently 
required. GPS and VLBI data from the DSN may be combined to generate a pre- 
cise UT1 series , while simultaneously reducing the time and effort the DSN must 
spend on platform-parameter calibrations. This combination is not straightforward , 
however , and a strategy for the optimal combination of these data is presented and 
evaluated. It appears that , with the aid of GPS, the frequency of required VLBI 
measurements of Earth orientation could drop from twice weekly to once per month. 
More stringent real-time Earth-orientation requirements possible in the future would 
demand significant improvements in both VLBI and GPS capabilities, however. 


I. Introduction 

One of the largest error sources for high-precision space- 
craft tracking and navigation is Earth-orientation 
variability. 1 ’ 2 The Deep Space Network (DSN) is cur- 
rently engaged in monitoring platform-parameter variabil- 


1 R. Treuhaft and L. Wood, “Revisions in the Differential VLBI 
Error Budget and Applications for Navigation in Future Missions,” 
JPL Interoffice Memorandum 335.4-601 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, December 31, 1986. 

2 S. W. Thurman, “DSN Baseline Coordinate and Station Location 
Errors Induced by Earth Orientation Errors,” JPL Engineering 
Memorandum 314-488 (internal document), Jet Propulsion Labo- 
ratory, Pasadena, California, July 25, 1990. 


ity through the Time and Earth Motion Precision Ob- 
servations (TEMPO) program of twice-weekly very long- 
baseline-interferometry (VLBI) measurements. This bur- 
den on the DSN radio telescopes is expected to increase as 
ever higher Earth-orientation accuracy is needed for space 
missions in the late 1990s and beyond. 

Earth-orientation measurements can be made using a 
Global Positioning System (GPS) receiver network cur- 
rently being installed at the three DSN sites and at a 
number of additional sites around the globe [1]. If an au- 
tomated, rapid-turnaround data collection and reduction 
system is in place, these GPS data can be used in conjunc- 


1 



( 1 ) 


tion with DSN VLBI data to produce a highly accurate 
Earth-orientation series in near-real time [2]. Such a sys- 
tem would help to reduce the frequency with which VLBI 
Earth-orientation measurements must currently be made, 
thereby allowing more time for direct spacecraft tracking 
and communications. 

The component of Earth orientation that is the most 
variable and unpredictable is the Earth’s rotation rate, as 
measured by Universal Time 1 (UT1). Its unpredictabil- 
ity is a major motivation for the twice-weekly TEMPO 
measurements currently being performed. GPS techniques 
promise to deliver daily estimates of the UT1 rate of 
change (also known as length of day or LOD) with a pre- 
cision of better than 0.1 msec/day [2]. Unfortunately, 
random error in LOD measurements can accumulate to 
produce estimates of UT1 (i.e., integrated LOD) that are 
in error by an amount that could exceed mission require- 
ments for UT1 accuracy. (Currently, a 50-nrad tracking 
accuracy requirement necessitates knowledge of Earth ro- 
tation at the ~0. 6-msec level; future missions may require 
0.1-msec- level real-time UT1 knowledge. 3 ) Hence, some 
sort of synergism is required between the GPS and VLBI 
techniques such that VLBI data would constrain long-term 
excursions in the GPS-derived UT1 time series while the 
GPS data would guarantee accurate monitoring of high- 
frequency UT1 variability. 

This article presents a method for combining DSN 
VLBI and GPS data. Section II discusses the nature of 
UT1 and LOD and highlights the need for frequent, ac- 
curate, LOD measurements. Section III reviews the ca- 
pabilities of both VLBI and GPS for measuring UT1 and 
LOD. Section IV discusses the combination of these two 
data types and various issues that need to be considered 
in order to generate the most effective combination. 

II. The Time Variability of UT1 and LOD 

The angle of rotation of the Earth is represented by a 
time, UT1 (equal to the rotation angle times a constant), 
and is usually discussed with respect to a reference rotation 
angle that is growing at a constant, predefined rate. Two 
commonly used reference angles, both based on atomic 
time standards, are Universal Coordinated Time (UTC), 
which is always maintained within 1 second of UT1, and 
International Atomic Time (TAI). The LOD is a measure 
of the rate of change of UT1, 4 and is defined by 

3 Treuhaft and Wood, op cit. 

4 LOD as defined in this article is more accurately termed A LOD, 
the excess length of day (the difference between the length of the 
day and Lob but the term LOD is retained based on common 
usage. 



where L is the excess length of day, U represents 
UT1-UTC, and L 0 is the nominal length of the day 
(86,400 sec) [3]. LOD currently has a magnitude of about 
2 msec, indicating that the value of UT1-UTC decreases 
by about 2 msec every 24 hours. 

LOD can vary in an unpredictable manner, typically 
changing by a few hundredths of a millisecond over one 
day. Figure 1 shows a sample of smoothed LOD over 
a recent 1-year period; note both the stochastic, high- 
frequency behavior and the more pronounced regular vari- 
ations at semi-annual and ^40-day periods. In this and 
subsequent figures illustrating LOD, tidal effects have been 
removed according to the model of Yoder et al. [4]. 

If unmodeled, this variability in LOD can, within a few 
days, lead to an accumulated UT1 error that exceeds the 
current accuracy requirement of 0.6 msec. To illustrate 
this, Fig. 2 shows the growth in error in UT1-UTC as a 
function of time for three simple cases for the behavior of 
LOD, according to 


U err or — Utrue U predicted — /. [~ r A(^’)]dr 


Utrue 


where r = time/Lo is dimensionless. In all three cases, the 
assumed reference value of UT1, U pre( n cie< i f is based on a 
value of LOD, L ri that remains constant with time; thus, 
U predicted varies linearly with time. In case one, LOD itself 
grows linearly with time, i.e, L(t) = L r (/ 0 ) + <*(f - f 0 )- In 
case two, LOD grows linearly with time for 2 days (as in 
case one), and then remains constant at its new value. 
In the third case, LOD is assumed to grow in propor- 
tion to the square root of time, i.e., L(t) = L r (to) -f cx(t 
— t 0 ) l/2 . The constant a is assumed to be 0.07 msec/day 
in the first two cases and 0.07 msec/>/day in the third 
case. (Note that in these examples, the units of time for 
the integration are given in days, while LOD and UT1 
are measured in milliseconds.) From Fig. 2, it is apparent 
that the maximum permissible UT1 error, 0.6 msec, -is ex- 
ceeded after 4 to 6 days. Although LOD behavior tends 
to be irregular, occasional secular variations in LOD (see 
Fig. 1) can result in unpredictable rapid growth of UT1 
error akin to that illustrated by the three cases shown in 
Fig. 2; at times such as these, frequent measurements of 
LOD or UT1 are critical. 
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A more realistic model for the day-to-day variability of 
LOD and UT1 may be found from power spectral den- 
sity (PSD) plots of LOD, such as that shown in Fig. 3. 
For periods longer than about 10 days, the spectrum 
is characterized by a (frequency)" 2 dependence, corre- 
sponding to a random- walk stochastic process [5]. LOD 
can thus be. characterized, over subseasonal time scales, 
as an integrated white-noise (i.e., a random-walk) process, 
while UT1 can be viewed as an integrated random walk. 
The white-noise PSD, Q, obtained from Fig. 3 is Q « 
0.005 msec 2 /day. After a time At = t ~ to, the expected 
variance of LOD is QAt, while the corresponding expected 
error is ~ 0.07\/A t msec (when At is measured in days). 
This stochastic model for LOD still produces an unaccept- 
ably large UTl error after 6 days. 

Given the typical variability of LOD, accurate UTl or 
LOD values must be available at least every 5 days in order 
to meet the current DSN requirements. In fact, the larger 
but less frequent LOD variations argue for measurements 
at least every 4 days. Furthermore, these time intervals do 
not take into account the required data processing time. 
The turnaround time for processing TEMPO data is typi- 
cally two to three days, at which time the UTl uncertainty 
has already grown to at least 0.2 msec. Since a new mea- 
surement must be available within 3 more days and it takes 
2 to 3 days. to. process the next data point, a new TEMPO 
measurement must be made within one more day. The 
net effect, assuming a 2- or 3-day turnaround time, is to 
require UTl measurements at least every 3 days. 


III. Measuring UTl and LOO With VLBI 
and GPS 

Very long-baseline interferometry is sensitive to the ori- 
entation of the VLBI baseline in a celestial reference frame 
defined by quasars. Since the angle corresponding to UTl 
is- one component of- the Earth’s orientation, in inertial 
space, UTl is a directly measurable quantity, although 
multiple baselines are required to unambiguously remove 
the effects of polar motion. Twice-weekly TEMPO mea- 
surements, one per week on each of two different intercon- 
tinental baselines, provide regular UTl and polar motion 
data to the DSN with typically a 2- or 3-day turnaround 
time (although 1-day turnaround has been demonstrated). 
Other VLBI networks, such as the National Geodetic Sur- 
vey’s International Radio Interferometric Surveying (IRIS) 
network and the U.S. Naval Observatory’s VLBI network 
(NAVNET), also provide regular UTl measurements, al- 
beit with slower turnaround times that are unsuitable for 
the needs of the DSN. The formal errors listed for the 
TEMPO data are currently at the 0.1- to 0.2-msec level 


[6], although combining and smoothing data from the two 
baselines can generate improved UTl estimates with for- 
mal errors below 0.1 msec. 

Each VLBI network relies on its own reference frames, 
both celestial and terrestrial, to produce Earth-rotation 
measurements. These reference frames may have small 
offsets with respect to one another, resulting in systematic 
differences between the Earth-rotation time series gener- 
ated by each network [7]. These UTl offsets are generally 
small and slowly varying, but they need to be considered 
in any precise estimate of UTl, especially if data from 
multiple techniques are combined. 

GPS methods differ fundamentally from VLBI tech- 
niques. To determine Earth orientation with GPS, one 
first needs to define a terrestrial reference frame. This 
is usually done by fixing the locations of a few select 
ground receivers. These sites are known as fiducial sites 
and are tied by local ground surveys to nearby, collocated 
VLBI antennas whose relative positions are known pre- 
cisely through VLBI [8]. A precise set of initial values, 
derived from VLBI and from satellite laser ranging (SLR), 
is used to orient this Earth-fixed frame with respect to 
a celestial reference frame, which then permits the GPS 
satellite orbit epoch states to be constrained in inertial 
space. The behavior of the satellite network is governed 
by dynamical models, such that movements of the solid 
Earth within that framework, i.e., Earth orientation, can 
be observed. UT1-UTC, however, cannot be determined 
simply by monitoring the relative motions of the terrestrial 
network and the GPS satellite orbits, as it is defined with 
respect to an arbitrary location on the celestial sphere. 
The signature in the data residuals of an error in UTl is 
identical to that due to an error in the satellites ascending 
node locations. In addition, since the GPS satellites have 
almost identical orbital radii, inclinations, and eccentrici- 
ties, the oblateness of the Earth induces precession of the 
GPS orbits that is nearly indistinguishable from a slow 
UTl variation. 

Although GPS techniques cannot directly measure the 
absolute value of UTl and have difficulty distinguishing 
long-term UTl changes from precession of the GPS orbits, 
they are sensitive to short-term changes in UTl. The large 
number of GPS satellites and receivers provide robust and 
redundant networks for measuring the displacements of the 
terrestrial frame with respect to the GPS orbits that result 
from rapid UTl variations. Regular GPS measurements of 
LOD are, in fact, already being made, although with a less 
than ideal satellite configuration and station network [9]. 
These values are weekly averages and have typical errors of 
about 0.15 msec. The capability of GPS to measure LOD 
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is expected to improve dramatically by the mid-1990s. The 
full constellation of 24 satellites will be available, along 
with a worldwide, high-precision GPS tracking network 
with the potential for near- real- time data processing. LOD 
measurements are thus expected at the 0.05-msec level or 
better in the next few years [2]. 

The terrestrial reference frame used in GPS processing, 
although tied to a VLBI terrestrial frame through the fidu- 
cial technique, is not necessarily identical to that of VLBI, 
since local site ties between VLBI and GPS antenna phase 
centers may contain errors. These site ties can also vary 
over time, due either to changes in antenna positions or 
to improvements in the local ground surveys. In addition, 
as VLBI site location and site velocity estimates improve, 
both VLBI and GPS reference frames and the Earth ori- 
entation time series may be affected. 

VLBI is a proven, currently operating technique to mea- 
sure UT1, whereas precise, routine, operational determi- 
nation of UT1 using GPS techniques lies at least a year or 
two in the future. Why, then, is there interest in supple- 
menting VLBI with GPS? 

As the tracking and navigation duties of the DSN ex- 
pand and higher precision is needed, more frequent UT1 
data will be required. As shown in Section II, the current 
frequency of TEMPO measurements is barely adequate 
in providing near-real-time Earth-rotation information; in 
fact, the TEMPO data must be supplemented with meteo- 
rologically derived atmospheric angular-momentum infor- 
mation (which has been shown empirically to be highly 
correlated with LOD) to ensure reasonable accuracy of 
current UT1 estimates [10]. The higher-precision future 
Earth-rotation requirements of the DSN will require more 
frequent VLBI data if alternative techniques are not avail- 
able. Dedicating DSN radio telescopes to obtaining Earth 
rotation at daily or 2-day intervals would heavily tax the 
already-overburdened DSN tracking schedule. Hence, fu- 
ture high-precision Earth-platform determination is crit- 
ically dependent on new high-accuracy technologies with 
the potential for rapid data turnaround, such as GPS. 


IV. Combining VLBI and GPS Data 

Three main issues need to be addressed in order to 
accurately and reliably combine GPS and VLBI Earth- 
rotation data. They are (1) growth of UT1 error with 
time, (2) re fere nee- frame definition and compatibility, and 
.(3) systematic technique-dependent errors. They will be 
discussed, in turn, below. 


A. Growth of Error in UT1 

Any strategy for the synergistic use of GPS and VLBI 
Earth- rotation data to monitor UT1 must take into ac- 
count error characteristics of both GPS and VLBI, along 
with the observed variability of Earth rotation. Since 
VLBI is able to measure UT1 directly whereas the strength 
of GPS lies in its ability to monitor LOD (the change from 
day to day in UT1), VLBI is needed to provide an initial es- 
timate of UT1. Subsequent measurements of LOD by GPS 
can then permit the calculation of UT1 by daily integra- 
tion. Since errors in UT1 will accumulate over time, VLBI 
data are also needed to calibrate the GPS-determined time 
series. This periodic VLBI calibration should overcome 
various problems related to long-term monitoring of Earth 
rotation with GPS, including the effects of orbit preces- 
sion, data noise, satellite force model errors, etc. 

Figure 4 schematically illustrates this plan for GPS and 
VLBI synergism. At some epoch, a VLBI observation is 
made, yielding a precise estimate (and formal error) of 
UT1. Daily measurements of LOD are then obtained by 
a rapid-turnaround global GPS network, and UT1 and its 
uncertainty determined by integrating these LOD values. 
At some point, the accumulating uncertainty in UT1 ne- 
cessitates a new VLBI observation, which fixes the value 
of UT1 at the new epoch and constrains its uncertainty. 
Note that Fig. 4 illustrates only a near- real- time method 
for estimating UT1; a UT1 time series of uniformly high 
precision will still be available after a delay of 2 weeks or 
more. 

To determine the frequency with which periodic VLBI 
measurements should be made, the sources of UT1 error 
must be evaluated. Formal errors for UT1 derived from 
TEMPO measurements are currently at the 0.1-msec level, 
although VLBI system improvements should reduce this 
uncertainty in the next few years. In this study, a formal 
VLBI-derived UT1 error of either 0.1 msec (conservative) 
or 0.05 msec (optimistic) and a 2-day turnaround time for 
processing these data are assumed. A UT1 measurement 
is assumed to represent the average value of UT1 over a 
few hours, obtained from back-to-back, 3-hour TEMPO 
VLBI sessions on each of two baselines. (Note that this 
differs from the current observing strategy of one 3-hour 
session on alternate baselines every few days.) 

Based on extensive covariance analyses [2], the LOD 
error anticipated from GPS when an optimal satellite and 
station configuration is established is less than 0.05 msec. 
To be conservative in this analysis, and for ease of presen- 
tation, we assume GPS-derived LOD also to be accurate 
to either 0.1 msec or 0.05 msec. Note that the GPS observ- 
able is, in actuality, an estimate of the UT1 rate-of-change 
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over a 12- to 24-hour period, but this is essentially equiv- 
alent to the average LOD. GPS data are assumed to be 
generated either every 24 or 48 hours, and to be available 
after a 24-hour processing delay. 

All error estimates mentioned are the one-sigma formal 
errors and are assumed to be Gaussian. In addition, no 
correlation is assumed between VLBI and GPS errors, even 
though technically there may be some correlation if any 
VLBI-derived parameters are used to initialize the GPS 
data processing scheme. 

Figure 5 illustrates a few examples of UT1 error be- 
havior. (The Appendix describes the mathematical model 
used to generate this and the subsequent figure.) Both 
the VLBI UT1 and daily GPS LOD measurements are 
assumed to have identical formal errors (either 0.10 or 
0.05 msec). At t = 0 days, a VLBI measurement of 
UT1 is made. The reduced data point is available at t = 
2 days; by this time, however, the real-time knowledge of 
UT1 based solely on VLBI will have degraded due to the 
stochastic nature of LOD. Fortunately, GPS data are be- 
ing obtained each day, providing a realistic assessment on 
day 2 of the LOD for day 0 and day 1 (recall that there is a 
24-hour delay for GPS data). Subsequently, GPS-derived 
LOD data with their formal errors control the growth of 
uncertainty in UT1. At some point (in this case, t = 
40 days), a new VLBI measurement is made in response 
to the accumulating uncertainty in UT1. Again, it takes 
2 days for the measurement to become available, and the 
cycle continues. 

In the case where formal errors for VLBI UT1 and GPS 
LOD are both 0.1 msec, the 0.6-msec DSN threshold for 
UT1 is reached in 35 days. Thus, DSN VLBI platform- 
parameter measurements would be needed approximately 
monthly. If the measurement errors for the two techniques 
drop to the 0.05- msec level, VLBI measurements of UT1 
might be needed only once every 3 or 4 months. In this 
latter case, the DSN could even collect GPS LOD data ev- 
ery other day and still free up the VLBI system for almost 
2 months. 

If, however, DSN requirements for near-real- time Earth 
rotation drop to. the 0.2-msec level, the picture is not 
so pleasant. This requirement could not be met with 
0.1-msec quality data. Only with high-quality 
(^0. 05-msec) data from both VLBI and GPS, and GPS 
LOD measurements obtained daily, could a 0.2-msec re- 
quirement be met without VLBI calibrations for any ap- 
preciable length of time. In this case, DSN VLBI mea- 
surements would be needed about every 2 weeks. This 
is still less than half as often as TEMPO currently oper- 


ates, and the alternative, TEMPO-only option would re- 
quire TEMPO measurements to be taken daily. To meet a 
0.1-msec real-time UT1 requirement, even VLBI and GPS 
data with errors at the 0.05-msec level are inadequate. 
Measurement errors must be reduced to the 0.02-msec level 
or better, and DSN VLBI calibration measurements must 
be made at least weekly to meet this more stringent Earth- 
rotation requirement. 

Figure 6 addresses a subtle issue concerning the near- 
real-time need for UT1. On any given day, the estimate 
of UT1 for the preceding day is an empirical value based 
on combined VLBI and GPS data. To generate an actual 
real-time estimate of UT1, this value for the preceding 
day must be adjusted for the behavior expected of LOD 
over the additional day, and the error estimate corrected 
accordingly. Two cases are illustrated for this theoretical 
behavior: a random walk and worst case linear growth. 
The random-walk model is simply the stochastic behav- 
ior of LOD described in Section II, while the worst case 
scenario corresponds to an unusually large daily change in 
LOD. These two cases correspond to growth in the uncer- 
tainty of LOD over one day of 0.07 msec and 0.10 msec, 
respectively. 

Figure 6 shows the expected error for these two mod- 
els if GPS data are taken either daily or, for purposes of 
comparison, every other day. A VLBI UT1 measurement 
with a 0.05-msec formal error is made on day 0 and be- 
comes available on day 2. GPS LOD measurements with 
a 0.05-msec formal error are made either daily or on al- 
ternate, even-numbered days. After every LOD measure- 
ment, the uncertainty in LOD is assumed to grow until the 
next LOD measurement, according to one of three models: 
no growth, random- walk (stochastic) growth, or rapid lin- 
ear (“maximum”) growth. The alternate-day LOD data 
curves oscillate because the uncertainty in UT1 24 hours 
after a GPS measurement is made is smaller than the UT1 
uncertainty 48 hours after the preceding measurement is 
made. Shown for comparison are the estimated uncertain- 
ties in UT1 if no GPS data were available to constrain the 
variability of LOD (the “VLBI only” curves). 

When daily GPS LOD data are used to estimate the 
current value of UT1, the modeled behavior of the LOD 
does not significantly affect real-time UT1 estimates on 
the day after the last data are taken. If GPS data are ob- 
tained on alternate days, however, the estimated UT1 un- 
certainties are significantly affected by the assumed LOD 
behavior. In this case, rapid LOD variations can produce 
UT1 errors exceeding 0.2 msec very quickly. Hence, once- 
a-day estimates of LOD by GPS are still the preferred 
measurement frequency. 
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B. Reference-Frame Compatibility 

As discussed in Section III, GPS and VLBI may gener- 
ate estimates of Earth-rotation parameters using slightly 
different terrestrial reference frames. Each frame should 
be internally stable over time, because a principal source 
of reference-frame instability, global tectonic motions, is 
a slowly varying process. Hence, differences between the 
two reference frames should manifest themselves in the 
Earth-rotation data as a more-or-less constant bias and, 
perhaps, a rate (i.e., a secularly growing offset) between 
the two time series. These bias and rate terms may be 
estimated empirically, based on previously measured time 
series of UT1 and LOD from the two techniques. The ter- 
restrial frames may also shift relative to each other in an 
unpredictable manner when fiducial sites are changed, an- 
tennas are moved, or improved local site ties between the 
GPS and VLBI antenna phase centers are determined. 

As long as the site network and data reduction method 
for each technique remain the same, the biases and rates 
should not significantly vary. However, changes in the 
reference site network, collocated site ties, VLBI radio 
source catalog, or data processing software may generate 
UT1 time-series glitches that need to be closely monitored. 
These time-series offsets or slope changes might be evident 
only weeks or months after they commence, when a suf- 
ficient quantity of data of diverse types has accumulated 
to distinguish the time-series characteristics of individual 
techniques. Fortunately, reference-frame adjustments are 
not expected to occur often. 

C. Systematic Errors 

Related to reference-frame discrepancies are the effects 
of systematic errors in each technique, either in data col- 
lection or in data processing, that contaminate the result- 
ing UT1 and LOD time series. Unlike reference-frame er- 
rors, these need not be of a long period. Systematic GPS- 
derived LOD errors have yet to be investigated, since the 
data themselves are not yet available at nearly the required 
frequency or precision. However, numerous orbit and base- 
line studies have isolated a number of possible systematic 
error sources for GPS techniques in general. These include 
orbit mismodeling, incorrect fiducial site ground ties, and 
mismodeled atmospheric propagation effects [11]. 

Isolating these systematic errors will require periodic 
campaigns of simultaneous GPS and VLBI Earth-rotation 
measurements. 5 Intercomparing these time series and in- 


5 The first of these multiple- technique simultaneous-observing ses- 

sions will occur in 199i. 


corporating the results of other Earth orientation mea- 
surement techniques, such as IRIS and NAVNET VLBI, 
SLR, and lunar laser ranging (LLR), should help to pin- 
point systematic errors. This type of analysis is ongoing 
for currently existing techniques [12]. 

V. Conclusions 

Within the next few years, high-precision GPS tech- 
niques will be able to generate Earth-rotation values of 
great utility to the DSN. By combining these LOD data 
with VLBI-derived estimates of UT1, a high-quality, real- 
time UT1 time series could become available for critical 
DSN spacecraft navigation and tracking applications. By 
reducing the amount of antenna time the DSN must com- 
mit to Earth-platform-parameter estimation, use of GPS 
data would permit a greater percentage of radio telescope 
time to be available for spacecraft telemetry and tracking. 
This would be all the more valuable given the expected 
increase in demand on the DSN radio antennas from the 
large number of planned and ongoing interplanetary mis- 
sions over the next decade. 

Although a high-quality GPS receiver network will al- 
ready be in place at the DSN sites and at numerous 
other locations around the globe in support of TOPEX/ 
POSEIDON tracking, the DSN must also establish a mech- 
anism for rapid GPS data collection and processing. A 
prototype system already exists for sub-daily orbit deter- 
mination, 6 but a fully operational system with 12-hour or 
less turnaround time and an explicit goal of measuring 
Earth orientation has yet to be implemented. 

Estimates of predicted UT1 uncertainty suggest that, 
given a small but reasonable GPS measurement error, the 
frequency of TEMPO VLBI measurements can be reduced 
from twice weekly to monthly or less (although this once- 
a-month TEMPO session must include measurements on 
both baselines). If a future real-time Earth-rotation re- 
quirement drops to a third of its present value (i.e., to 
<0.2 msec), the only feasible way to fulfill this require- 
ment, short of daily TEMPO measurements, is the syner- 
gistic combination of daily GPS LOD measurements and 
weekly or semimonthly VLBI UT1 measurements. Fur- 
thermore, to achieve a 0.1-msec real-time Earth-rotation 
capability, significant improvements in both VLBI and 
GPS techniques are also needed. Thus, both VLBI and 
GPS appear to be techniques essential for meeting DSN 
platform-parameter requirements in the coming years. 


6 G. Blewitt, personal communication, Jet Propulsion Laboratory, 
Pasadena, California, 1990. 
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Fig. 5. Estimated one-sigma error In real-time UT1-UTC as a function of time for three different 
measurement scenarios. Lines are used to connect the points only; they do not represent the 
estimated errors at intermediate times. Dotted lines indicate retroactive uncertainty in UT1 
following the VLBI measurement. 
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Fig. 6. Expected real-time one-sigma error in UT1-UTC when an LOD behavior model is in- 
cluded. Note that lines are used to connect the points only; they do not represent the estimated 
errors at intermediate times. 
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Appendix 

Modeling the UT1 Error Structure 


The mathematical model for the evolution of UT1 error estimates shown in Figs. 5 and 6 is presented below. 
Measured data are assumed to consist of the following: U m |o (UT1-UTC at time to) with formal error <Tu\ 0} and LOD 
estimates (L m \i) at times £{,1 = 0,...,/, with formal errors The true measurement errors in UT1 and LOD are 

denoted by cy | 0 and respectively. All times are normalized, i.e., £,* = time(i)/ L q. Recall that LOD estimates are 
actually estimates of —AUTl/At over a time interval At. 


The predicted value of XJT1 at t n , where t n > </, is 

i 

U P (tn) “ t^m|0 “ ^ ^ (^» ““ — 1 ) “ 

i= 1 

where L'(r) is the predicted behavior of LOD at normalized time r. In the absence of data, LOD is usually assumed 
to remain constant, i.e., L'(r) =£ m p for r > £/. The true value of UT1 is represented by 


/ 

Jti 




(A-l) 


U(t n ) = U 0 - f n L(r)dr = [t / m , 0 + £u|o] - /*' L(r)dr + f" L(r)dr 

J to J t o J t X 


(A-2) 


where Uo and L(r) denote the true UT1 and LOD values, respectively. 


The true value of LOD may be modeled within the above equation as 


/ L{r)dr = ^ f L(r)dr = Z |,_i(t 
Jto i-i i=1 


i-ti-i) 


(A-3) 


where L |i_i is the mean value of LOD over the interval £,*_i to t{\ it is assumed that L = L m + £L\i-\- The 
final term in Eq. (A-2), the effect of LOD variability following the last LOD measurement at time £/, is modeled as 


f L(r)dr= f [L'(t) + rj(r)]dt = / [X m( , + rj(r)]dr 

Jt t Jt X Jtj 


where the difference between the true and assumed LOD is 


(A-4) 


V( T ) = £l\i + 



(A-5) 


In Eq. (A-5), the first term on the right is the LOD measurement error at t [ , while the final term in this expression 
emerges from the random-walk nature of LOD. The variable w(/i) describes a white-noise stochastic process whose 
integral over time generates the random-walk variations of LOD. Alternatively, the true LOD may be modeled as an 
analytic function of time, i.e., rj(r) = d\i + f(r — ti). 
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The error terms and white-noise behavior in the above equations may be described by the following set of expectation 
values: 


(£t/|o £t/jo) = {<?u\ o) 2 

( e L\i ZL\j) = {?L\i) 2 Sij ► 

(u(ti)u(p)) = Q 6 (fi-p) J 


(A-6) 


while all other expectation values of u>(/i), £u\o, £l\{ and their products are assumed to be zero. In Eq. (A-6), 6{j is the 
Kronecker delta, 6(p — p) is the Dirac delta function, and Q is the PSD of the white-noise process. 

The error in UT1 as a function of time is given by 


/ 

£u(tn) = Up(t n ) - U(tn) = -£t/|0 + ^ [g£,|t-l(*» ~ *»-l)] 

t = l 

The variance of this error is 


rtn 

+ / V(r)dr 
Jti 


Var(ec/(^n)) = {£u{U)eu{U)) 


= (o-(/|o ) 2 + - ti-i ) 2 + [ f {v(r) v(p))dp dr 

{-1 J T — t\ Jp-U 


(A-7) 


(A-8) 


where most of the expectation values have been evaluated. The rightmost term in Eq. (A-8) consists of two parts, one 
corresponding to the error of the final LOD measurement, £l\i> and one corresponding to either the random-walk or 
analytic model for LOD. Thus 

Q(tn ~ f /) 3 

3 2 (A-9) 

[//,” f( T ~ 

The upper quantity on the right corresponds to the random-walk generated uncertainty, and the lower to the uncertainty 
due to the modeled analytic behavior of LOD. Figure 5 is generated using the first two terms on the right-hand side of 
Eq. (A-8), while Fig. 6 utilizes all of Eqs. (A-8) and (A-9). 


/• / 

J T — tl J P = tl 


( 7 ?( r ) n(p)) d P dT = (°T|/) 2 (fn - u) 2 4* < 
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